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1 Introduction 

CP violation plays an important role in particle physics since its discovery 40 years 
ago pQ . For more than 20 years this phenomenon appeared as confined in a particular 
sector of particle physics, through the mixing between states of opposite CP in the 
neutral kaons. The unambiguous discovery in the late 1990s, after the early indication 
by NA31 [2], of direct CP violation in K -> 2vr decay by the NA48 [3] and KTEV @] 
experiments and the discovery of CP violation, in its various forms, in the neutral 
B meson system [5] represented important breakthrough in the understanding of 
the particles dynamics. A complete as possible study of the tiny effects due to the 
violation of this symmetry in all the systems where it can be carried out, represents an 
important window on the contribution of new physics beyond the Standard Model: 
in fact new effects could appear, in particular in the heavy quark loops which are 
the core of the mechanism allowing the CP violation in the mesons decay. In the 
kaon sector the most promising places, besides e /e, where this kind of contributions 
could play some role are the rates of GIM suppressed rare decays and the charge 
asymmetry between charged kaons. In particular the K — > 3ir asymmetry could give 
a strong qualitative indication of the validity of the CKM description of the direct 
CP violation or reveal the existence of possible sources outside this paradigm. 

In principle Direct CP violation in K can be detected comparing the different 
decay amplitudes in the 3n mode: 

\A(K + -> tt+tt+tt-)! ^ \A{K~ -> 7^+01 

Experimentally the easiest way to study the difference between the charge conjugate 
modes is to compare the shape of the Dalitz Plot distribution instead of the decay 
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Asymmetry 


# of events 


Experiment 


Ag = (19 ± 125) ■ 10" 4 
A° q = (2 ± 19) • 10" 4 


115K 
620K 


CERN PS(1975) © 
Protvino IHEP (2005) [7] 


A c g = (-70 ± 53) ■ 10" 4 
A° q = (22 ± 15 ± 37) • 10~ 4 


3.2M 
54M 


BNL AGS (1970) [8] 
HyperCP (2000) prelim [9] 



Table 1: Summary of the experimental situation both in "neutral" (A®) and in 
"charged" (A c ) mode, before the NA48/2 results 

rates. The small phase space in the three pion decay mode allows to expand the 
matrix element in terms of the, so called, Dalitz variables u and v: 

_ _(s 3 - s ) (s 2 - Si) 



ml ml 

7T IT 

These variables are defined using the Lorentz invariant Si = (px — Pi) 2 , where px 
is the kaon four momentum and Pi are the pion four momenta [i = 1,2,3 the latter 
being "odd" pion) and s = (si + s 2 + s 3 )/3. Exploiting the Dalitz variables the 
matrix element can be written as: 

\M(u,v)\ 2 ~ 1 +gu + hu 2 + kv 2 + ... (1) 

where g,h,k are the linear and quadratic slope parameters. Using the fact that |g| >> 
\h\, \ k\, the CP violation parameter 

is defined using only the linear slopes given in ([T]). Being g + relative to K + decay 
while g~ to K~ , the parameter defined above is different from only if an asymmetry 
exists between the matrix element describing kaon decay of opposite charge. Theo- 
retical predictions for the A g parameter both in the K — ► ix^ti^ix {A™, the so called 
"neutral" mode) and in the K — > 7r ± 7r + 7r _ (Ag, the so called "charged" mode), are 
very difficult and the available predictions varying from 10~ 6 to few 10 -5 are unre- 
liable; calculations in the framework of theories beyond the standard model predict 
a substantial enhancement of this parameter up to the level of few 10~ 4 . Several ex- 
periments in the past have searched for asymmetry both in "charged" and "neutral" 
mode. The sensitivity reached so far is at level of 10~ 3 , as summarized in table [TJ 
The main goals of the NA48/2 experiment are to reach the sensitivity of 10~ 4 both 
in "neutral" and in "charged" mode, to investigate the possibility of non standard 
model contributions to the CP violation in charged kaon decays, thus covering the 
gap existing between the experimental results and the SM theoretical predictions. 
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Figure 1: The NA48/2 beam line and detector. Not in scale. 



The large amount of K — ► 37r collected by NA48/2 experiment, allows a very precise 
measurement of the Dalitz plot parameters and, as recently shown [10], the study of 
the neutral Dalitz plot density allows to extract important informations about the 
pion scattering length. 



2 Beams and detectors 

The simultaneous collection of K + and K~~ decay with the same apparatus is an 
essential point of the asymmetry measurement. A novel high intensity beam line was 
designed in the high energy hall (ECN3) at CERN SPS (fig. [Tj). The charged particles 
(kaons, pions, muons and electrons) are produced by 400 GeV high intensity protons 
beam(~ 7 ■ 10 11 protons per pulse), from the SPS accelerator, with a 40 cm long 
and 2 mm in diameter beryllium target with an angle of zero degrees. A magnetic 
device, called first achromat, selects the momentum of the beam in the range (60 ± 
3)GeV/c, splitting the two charges. After being recombined the beams are focused 
by a quadruplet of quadruples, before another splitting in the second achromat. The 
second achromat houses the first detector along the beam line, called KABES 
a spectrometer measuring the particles' momentum with a resolution of ~ 1% (this 
detector is used studying rare decays). The two charged beams, recombined again 
along the beam axis, contain ~ 6.4 • 10 7 particles per 4.8 s spill, with a charge ratio 
K + / K~ ~ 1.8 (irrelevant for the charge asymmetry measurement) and 12 times more 
pions with respect to the kaons: however the pions decay products remain into the 
beam pipe, because of the small momentum and do not cross the detectors. The 
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decay region is housed in an evacuated tube ~ 114m long and ~ 2m in diameter. 
The beams in the decay region are superimposed within 1 mm with a total width of 
~ 5mm, so that both K + and K~ decays illuminate with the same acceptance the 
same detector. The central detector is based on the old NA48 detector, described 
elsewhere p2]. For the asymmetry and Dalitz plot parameters measurement the 
two main detector are the spectrometer and the LKr calorimeter. The magnetic 
spectrometer works with a px = 120MeV/c kick and the resolution in momentum 
(GeV/c) is <J p /p = 1.0% + 0.044%. In order to manage the higher intensity with 
respect to the previous NA48 runs, the drift chambers read out has been rebuilt. In 
order to collect the gammas coming from the neutral pions decay a Liquid Krypton 
(LKr) calorimeter with a resolution in energy (GeV) of 

" 3 - 2% + ?* + 0.42% 



E E 

is employed. The very good resolution in the reconstructed kaon mass (1.7 MeV/c 2 
for K — > /T ± 7r + 7r _ and 0.9 MeV/c 2 for K tt ± it it ) allows a precise calibra- 
tion and monitoring of the characteristics and performances of the apparatus. The 
data collection is based on a multilevel trigger system. The first level (LI) uses the 
information coming from a plastic scintillator hodoscope and from a dedicated LKr 
readout, in which the number of peaks in the calorimeter's energy deposit are com- 
puted. The second level (L2) is based on processors for a fast DCH reconstruction. 
In particular the number of reconstructed vertexes with 2 or 3 tracks, is used to col- 
lect — > 7r ± 7r + 7r^ events and the missing one-track mass, assuming the nominal 
kaon momentum (60 GeV) and direction (z axis), is used to collect K — > 7r ± 7r 7r°, 
rejecting the main — > tt^tt background. NA48/2 collected data during two runs 
in 2003 and 2004, for a total of ~ 100 days of data taking. About 18 • 10 9 triggers 
have been registered on tape, for a total of more than 200 TB. 



3 Charge asymmetry measurement strategy 

The asymmetry method is based on the comparison between the u projection of 
the Dalitz plot distribution, in order to extract the difference between the matrix 
element linear components. This difference, defined as Ag = g + —g~, can be extracted 
considering the ratio between the density in the u distribution for K + and K~ decays. 
The ratio between the two distribution can be written as: 

N K+ 1 + (g + Ag)u + hu 2 

R{u) = —— oc — — — 3 

N K - 1 + gu + hu 

where g = (g + + g~)/2. From Ag the asymmetry parameter A g can be easily eval- 
uated using the relation A g ~ Ag/2g. The simultaneous collection of decays coming 
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from beams with similar momentum spectrum and a similar detection efficiency and 
acceptance of the decay products, are fundamental points to control the instrumental 
charge asymmetry. However, the presence of magnetic fields both in the beam sector 
(achromat) and in the detector (spectrometer magnet) could introduce an intrinsic 
charge dependent acceptance of the apparatus. To equalize this asymmetry the main 
magnetic fields were frequently reversed during the data taking. During the 2003 run 
the magnet spectrometer polarity was reversed every day while the achromat magnets 
polarities every week. In the 2004 run the reversal was more frequent: every about 3 
hours for the analyzing magnet and 1 day for the beam transport line magnets. It is 
possible to redefine the ratio in ([3]) taking into account the magnetic field alternation. 
For instance, for a given achromat polarities, the ratios: 

]V B ~ jv B+ 
Mu) = jfi , Rs(u) = -ft (4) 

are defined using the same side of the spectrometer, in the sense that the numerator 
and the denominator in these ratios contain particles deflected in the same direction. 
The subscripts J and S represent the particles bending direction according to the 
geographic position of the Jura (J) and Saleve (S) mountains, respectively on the 
left and right side, with respect to the NA48/2 beam line direction. Considering the 
possible achromat polarity, four independent ratios can be built exploiting the four 
different field combinations: instead of the single ratio ([3]), a quadruple ratio can be 
defined as: 

-f l + (g + Ag)u + hu 2 \ 4 

R{u) = RusRujRdsRdj ~ R — —7-5 (5) 

y 1 + gu + nu z J 

where U and D stand for the path, up or down, followed by K + in the achromat 
system, and R is an inessential normalization constant. This method is independent 
on the relative size of the four samples collected with different fields configuration and 
on the K + and K~ flux difference. In the value of A g extracted from the quadruple 
ratio (jSJ), the benefits due to the polarity reversal are fully exploited and the main 
systematic biases due to instrumental asymmetries cancel out. In particular in the 
quadruple ratio there is a three fold cancellation: 

• local detector bias (left-right asymmetry), thanks to the fact that each single 
ratio is defined in the same side of the detector; 

• beam local biases, because in each single ratio the path followed by the particle 
through the achromat is the same; 

• global time variation, because the decays from both charges are collected at the 
same time (this is not true for the same single ratio in which the numerator and 
the denominator are collected in different period). 
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The result remains sensitive only to the time variation of the detector with a char- 
acteristic time smaller than the inversion period of the magnetic field, if this effect 
is charge asymmetric and u dependent. Other systematics biases induced by effects 
not canceled by the magnetic field alternation of the magnetic fields (for instance the 
Earth's magnetic field in the decay region and any misalignment of the spectrometer) 
have been carefully corrected. The intrinsic cancellation of instrumental asymmetries 
in the quadruple ratio allows to avoid the use of a MonteCarlo simulation. Neverthe- 
less a GEANT3 based MonteCarlo, including the full geometry description and time 
variations of beam characteristics, DCH inefficiency and spectrometer alignment, is 
used for systematics studies and as a cross-check of the result. 



The reconstruction in the "neutral" mode is based on LKr to construct the u variable 
and on the spectrometer to define the event charge and close the kinematics. The 
possibility to define the u variable using the 7r°'s only is a strong point in this kind 
of analysis, because of the charge independence of them. Anyway an alternative u 
reconstruction can be done using the DCH and the KABES informations, to obtain 
a result, useful as cross-check, with different systematic effects. The fiducial region 
of the detectors is chosen to avoid edge effects. In particular in order to symmetrize 
the small difference in the spectrometer acceptance between decays coming from K + 
and K~ beams, the spectrometer's inner radial cut is chosen according to the actual 
beam position, periodically measured as the average of the reconstructed transverse 
vertex position using three charged pion events. The decay vertex is reconstructed 
from the 7 impact point position on the LKr, for each 7T°, by using the formula: 



Among all the possible 7 pair, the correct pairing is chosen minimizing the difference 
between the two 7r°'s vertexes. The final vertex is obtained as arithmetic average. 
The kaon invariant mass is obtained including the charged pion measurement from 
the DCH, in order to reduce the events background by requiring \M K ± — Mpdg\ < 
6MeV/c 2 . A total of 59.3 • 10 6 K + and 32 • 10 6 K~ events has been selected; In 
plot [2] the reconstructed Dalitz Plot is shown. The photon position on the LKr is 
corrected to take in to account the calorimeter projectivity. The measurement of 
the charged momentum, slightly biased by variable DCH misalignment, is corrected 
exploiting the condition M K + = M K - in the K — > 7i ± 7i + 7i^ . Using the same decay 
mode the magnetic field inversion is monitored online at level of 10 -3 studying the 
reconstructed kaon mass with respect to the nominal (PDG) kaon mass. The effect 
of the residual magnetic field in the decay region (the so called Blue Field), mostly 
due to the earth magnetic field, is corrected using the maps obtained from a direct 



3.1 Results in K 





(6) 
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Figure 2: Dalitz plot distribution for 



Figure 3: Dalitz plot distribution for 

— > 7r ± 7r + 7r _ 



measurement before the data taking. Several sources of potential systematic bias 
have been considered. The effect due to the acceptance has been evaluated studying 
the stability of the result varying the cuts definition. The MC has been used to 
study the contribution of the pion decay in flight to the total systematics error. 
The contribution to the systematics of the online trigger system has been carefully 
studied. The efficiencies of the LI and L2, have been evaluated using data collected 
with control triggers uncorrelated with the main trigger. In the neutral mode the 
LI is essentially obtained with a coincidence between a signal (Ql) coming from 
the scintillating hodoscope and a signal (NTPEAK) compatible with a deposition of 
four clusters in the LKr. The L2 (MBX) is based on the algorithm that rejects the 
tt^it . The main source of systematics comes from the neutral part of the LI trigger 
(NTPEAK) being limited by the number of events in the control sample. In table [2] 
a summary of systematics is shown. The preliminary result in the slope difference 
using the whole statistics is: 

Ag = (2.7 ± 2.0 stat ± 1.2 syst ± 0.3 ext ) x 1CT 4 

where the external error is due to the ~ 3% error [13] on the knowledg^j] of the g 
value. The resulting charge asymmetry parameter is: 

A™ = (2.1 ± 1.6 stat ± l.0 syst ± 0.2 ext ) x lCr 4 = (2.1 ± 1.9) x 10" 4 

This result is fully compatible with the SM prediction and is almost one order of 
magnitude better than the previous measurements [B] [7]. 

lr rhis contribution becomes negligible using the new go measurement |14j 
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±0.1 
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±0.8 


L2 Trigger 


±0.6 


Total systematic uncertainty 


±1.2 



Table 2: Systematics in K — * 7r ± 7r°7r° 
3.2 Results in A' ± ^tt+tt 

The offline reconstruction of the K ± — > tt ± 7c + tt~ is totally based on the Spectrometer. 
The decay vertex is obtained extrapolating the track segment from the first two 
chambers from the spectrometer to the decay volume, taking into account the presence 
of the Blue Field in the decay region. The track momentum is rescaled, to compensate 
the variation of the DCH alignment, as described above. In the three charged pion 
mode the chambers' acceptance and the spectrometer performance are most critical 
with respect to the neutral case where the spectrometer is used only to tag the event. 
In particular the beam pipe crossing the chambers determines the main difference in 
the acceptance between the two beams for which the beam optic can not control the 
transverse position better than ~ 1mm. The cut centered on the actual beam position 
(at level of the first and last chamber) must be applied to all the pions, resulting in 
a reduction of ~ 12% of the whole statistics. In plot [3] the reconstructed Dalitz Plot 
is shown. The only relevant physical background come from the pion decay in flight. 
More than 2 • 10 9 K + and 1.1 • 10 9 K~ decays are collected in this channel. The main 
systematics in the charged mode comes from the pion decay in flight as reported in 
table [3] where the contributions to the systematic error are summarized. A simpler 
fitting function can be used in the charged mode case with respect to (JSJ) to extract 
the Ag value, exploiting the relative smallness of the measured g: 

R(u) oc 72(1 ± Aguf 

The preliminary result, based on the 2003±2004 data taking, is: 

Ag = (0.6 ± 0.7 stet ± 0.7 syst ) x 10~ 4 
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Table 3: Systematics in K ± — > 7r ± 7r + 7r 

leading to a charge asymmetry parameter of 

A c g = (-1.3 ± 1.5 stoi ± 1.7 syst ) x 10- 4 = (-1.3 ± 2.3) x 10" 4 

Also in this case the goal to increase by a factor 10 the sensitivity with respect to 
the previous measurement has been reached. The reason for a similar precision of 
asymmetry results in "neutral" and "charged" mode, in spite of the different statistics, 
lies in the fact that the Dalitz Plot density is most favourable in the "neutral" mode. 
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4 Dalitz plot parameters measurement in K 

In fig. H] the 7r°7r° invariant mass square (Mq ) distribution (proportional to the u 
distribution) is shown. The change of slopes, seen for the first time by NA48/2 at 
Mq = (2m n +) 2 , can not be explained by the simple matrix element parametrization 
given in ([T]). This structure has been interpreted by Cabibbo [TU] as due to the 
rescattering process 7r + 7r~ — > 7r°7r° coming from the K ± — > n ± n + n^ decay. In fact the 
K ± — > tt^tc ^ amplitude can be written as the sum of two terms (just considering the 
first rescattering order): the direct 7r°7r° emission A4q, parametrized by the standard 
polynomial expansion, and the terms due to the rescattering process This last 
term, proportional to the difference (ao — (12) between the pionic scattering length for 
1=0 and 1=2, is real below the threshold of 2m n + and imaginary above. The total 
amplitude can be written as: 

U 2 _ J (-M0) 2 + (A-li) 2 + 2 ■ A-loAfi s re < Aml + 

M I (M ) 2 + \Ml\ s^>Aml + [() 
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Figure 4: U distribution for K — > t^tt :^ . The 7r°7r° invariant mass square (Mq ) is 
related to the u variable by u = (Mq — s )/m^. In the arrow position (corresponding 
to (2m 7T +) 2 ) it is possible to see the "cusp" structure. 



The AioAii term gives a destructive interference below threshold. Other rescattering 
diagrams can be included in a systematic way as shown by Cabibbo and Isidori |15j . 
Fitting data with only the M., term yields a fair agreement only above the 2m n + 
threshold, because of the anomaly introduced by the pions strong rescattering in the 
Mq distribution, while the standard expansion is not enough to explain the complex 
dynamics contributing to the whole decay amplitude. Using the Cabibbo one-loop 
model, the fit quality increase giving a x 2 °f 420.1 for 148 degree of freedom, as 
shown in fig. [6] . Including the two-loop Cabibbo-Isidori approach the x 2 becomes 
more reasonable (158.8 for 146 degree of freedom). Near the threshold the relative 
pions velocity decreases and the possibility to have electromagnetic it + 7t~ bound 
states increases. However the description of the so called pionic atoms (pionium) 
needs particular care, because of the Coulomb interaction correction and the critical 
experimental resolution (the fit obtained including the pionium is shown in the third 
plot in fig. EJ). For that reason we prefer to exclude 7 bins around the threshold 
position to perform the final fit (last plot in figJH] ) in which we have a x 2 °f 145.5 
for 139 degree of freedom. The results [16] . based on 23- 10 6 events collected in 2003, 
are obtained setting k, the quadratic v slope, to 0: 

g° = 0.645 ± 0.00A stat ± 0.009 s ^ 

h' = -0.047 ±0.012 stet ±0.011 fl „ st 
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Figure 5: A = (data — fit) /data is 
shown function of M^o^o. The cusp 
threshold is indicated by an arrow. The 
fit x 2 is reasonable only if the fit is re- 
stricted to the indicated fit interval. 
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Figure 6: The application of the rescat- 
tering hypothesis improves the quality 
of the fit (see text) 





NA48/2 results 


PDG06 


g 
h 
k 


(-21.131 ± 0.009 stat ± 0M2 syst )% 
(1.829 ± 0.015 stat ± 0.036 SJ/St )% 
(-0.467 ±0.005 stat ±0.011 SJ/st )% 


(-21.54 ±0.35)% 

(1.2 ±0.8)% 
(-1.01 ±0.34)% 



Table 4: Results for Dalitz Plot parameter in charged mode. 



where ti + (1/A)g 2 = h. The data are compatible with a non zero value for the k 
parameter, never measured before. The value of the fit 

k = 0.0097 ± 0.0003 stat ± 0.0008 SJ/St 

is still preliminary. The (do — a 2 ) value is not affected by the k ^ term, but the g 
and h values change respectively by ~ 2% and ~ 25%. 

5 Dalitz plot parameters measurement in K ± — > 

Thanks to the huge statistics collected in K ± — > 7i ± 7i + n^ decay and to the well tuned 
MC, very precise measurement of the Dalitz plot parameters can be performed. The 
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pion rescattering effects can inffuence the matrix element also in the "charged" mode 
but, being on the border of the Dalitz Plot, is not so evident like in the "neutral" 
mode. The goal of the first and preliminary study presented here is to measure the 
parameters in the standard polynomial expansion verifying the validity of (CD). The 
parameters (g,h,k) are obtained minimizing the 



{Fdata — NF MC 



2/ 7 7 »n data 1,1 MC 

X (g, h, k, N) = , p2 , N2F 2 



2 

data iv 1 MC 



where F represents the population (in data or in M C) in the (u,v) bin. The MC pop- 
ulation is obtained adding the 4 components corresponding to the four possible terms 
in pi) . The relative weights, obtained from the fit, are the polynomial expansion 
parameters. The coulomb correction is applied to take into account the pion electro- 
magnetic interaction. The main contributions to the systematic uncertainty, at this 
analysis stage, come from the momentum scale in the charged pions measurement, 
the kaon momentum spectrum in the MC and trigger inefficiencies. The preliminary 
result is based on 0.47 ■ 10 9 events collected in the 2003 run. In table H] the results 
are presented and the agreement with the PDG06 values is shown. The previous 
measurements, by experiments made in 70s, are one order of magnitude less precise 
with respect to the NA48/2 measurement, based on 1/4 of the whole statistics. 



6 Conclusions 

The main goal of the NA48/2 experiment was to measure, with a precision at level 
of 10 -4 the charge asymmetry parameter A g , both in K — > 7T ± 7r + 7T _ and in K — > 
^ti ^ decays. The preliminary results obtained in the "neutral" (Ag) and "charged" 
(Ag) mode: 

A n g = (2.1 ± 1.6 stat ± 1.0 syst ± 0.2 ex t) x 10- 4 = (2.1 ± 1.9) x 10~ 4 

A° g = (-1.3 ± l.5 stat ± 1.7^) x HT 4 = (-1.3 ± 2.3) x KT 4 

are compatible with the SM predictions and with our previous results based on partial 
samples [T7] [T8] . 

In K — > 7r ± 7r°7r the standard polynomial matrix element expansion is not enough 
to describe the observed tt°tt° invariant mass spectrum. Taking into account the 
rescattering processes, whose contributions are proportional to the (ao — a 2 ) term the 
corresponding slope are found to be (setting k=0) 

g° = 0.645 ± 0.004 stat ± 0.009^ ti = -0.047 ± 0.012 stat ± 0.011^ 

However a ~ 1% for the k quadratic slope is obtained with a complete fit (preliminary 
result). 
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In the charged mode the standard polynomial fit has been employed. The Dalitz 
plot parameters have been remeasured with higher precision with respect to the pre- 
vious old measurement. 
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